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We examine the scattering properties of a slowly and periodically driven mesoscopic sample using 
the Floquet function approach. One might expect that at sufficiently low driving frequencies it 
is only the frozen scattering matrix which is important. The frozen scattering matrix reflects the 
properties of the sample at a given instant of time. Indeed many aspects of adiabatic scattering can 
be described in terms of the frozen scattering matrix. However, we demonstrate that the Floquet 
scattering matrix, to first order in the driving frequency, is determined by an additional matrix 
which reflects the fact that the scatterer is time-dependent. This low frequency irreducible part of 
the Floquet matrix has symmetry properties with respect to time and/or a magnetic field direction 
reversal opposite to that of the frozen scattering matrix. We investigate the quantum rectification 
properties of a pump which additionally is subject to an external dc voltage. We split the dc current 
flowing through the pump into several parts with well defined properties with respect to a magnetic 
field and/or an applied voltage inversion. 

PACS numbers: 72.10.-d, 73.23.-b, 73.40.Ei 



I. INTRODUCTION 

The interplay of quantum mechanical interference with 
quantized energy exchange results in a quantum pump ef- 
fect which is investigated intensively both experimentally— - — 
and theoreticallyi- - — This phenomenon being promising for 
manipulating and controlling the passage of electrons through 
mesoscopic circuits is of fundamental interest. Adiabatic driv- 
ing involves only low energy exchange and avoids excitations 
into inelastic channels which degrade the quantum proper- 
ties of the system. In this work we investigate the magnetic 
symmetry properties of the dc-current of a quantum pump 
which might operate in the presence of applied voltages and 
temperature gradients. 

The experimentally measured adiabatically pumped dc 
current 1 flowing through a chaotic cavity with periodically 
varying shape is symmetric in magnetic field H. That is 
in seeming contradiction with the theorjiSi— predict- 
ing that the pumped current has no definite symmetry under 
magnetic field reversal. As a result it was conjectured— - — 
that the current measured in Ref.Qis caused by a classical rec- 
tification effect. Indeed subsequent measurements 3 confirmed 
that for slow one-parameter driving there is a symmetric in 
magnetic field induced current whose origin is classical rec- 
tification. Nevertheless one can not exclude the possibility 
that the current measured in Ref. [l| contains also the con- 
tribution coming from the quantum pump effect. To check 
it, perhaps, it is necessary to investigate the system in a less 
symmetric setup, i.e., with reservoirs having different elec- 
trochemical potentials or temperatures. In the present paper 
we give a simple example when the pumped current has or 
has not an odd in magnetic field contribution depending on 
whether there is or there is no applied voltage. Further exper- 
imental and theoretical efforts to detect and distinguish the 



quantum pump effect are highly desirable in view of a possible 
application in quantum information processing devices—*— 

The aim of the present paper is to explore in detail the 
symmetry properties of the adiabatic current generated by 
the periodically driven mesoscopic conductor. To this end 
we represent the Floquet scattering matrix at low driving 
frequency uj as a sum of different terms with well defined 
symmetry properties (e.g., with respect to a magnetic field 
direction reversal). One term reflects the symmetry of a sta- 
tionary scattering process while the other term vanishing at 
uj — *• has symmetry properties opposite to a stationary scat- 
tering process. Based on such a representation we divide the 
dc current into parts with well defined symmetry properties. 
That opens up additional possibilities for the experimental 
detection of the quantum pump effect. 

In particular, in the two terminal case, we find a voltage 
dependent contribution to the pumped current which is odd 
in magnetic field. At small voltage this current is linear in V . 
Thus for small magnetic fields the dc-current has a component 
which is proportional to the product of frequency, magnetic 
field, and applied voltage. For comparison we recall that in 
the stationary case, for a two-terminal conductor, the current 
linear in voltage (or, alternatively, the conductance) is an even 
function of a magnetic field4»i— A current that is odd in mag- 
netic field appears only in the nonlinear voltage regime— *— 
and is caused by electron-electron interactions. In contrast, in 
the non- stationary case considered here even non-interacting 
electrons can show a response that is odd in magnetic field 
and linear in applied voltage. 

Recently the magnetic field symmetry of the dc current 
through an open quantum dot subject to a one-parameter 
potential oscillation has been investigated experimentally and 
theoretically as a function of frequency^ In contrast, in the 
present paper we consider a two-parameter oscillation and 
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investigate the magnetic field symmetry of the dc current in 
the presence of adiabatic parametric quantum pumping. 

The paper is organized as follows. In Sec[H]we briefly con- 
sider the Floquet function approach to scattering of electrons 
at a periodically driven mesoscopic conductor and analyze 
the consequences of microscopic reversibility. We introduce 
an exact representation for the scattering matrix at low driv- 
ing frequency ui. According to this representation the Floquet 
scattering matrix elements (up to linear in u) terms) are pro- 
portional to the elements of both the stationary scattering 
matrix So and a residual Floquet matrix A which exhibits 
symmetry properties opposite to those of So. The symme- 
try properties of So are dictated by micro-reversibility, and 
the residual Floquet matrix A reflects directly the breaking 
of these symmetries due to the driving of the sample. Using 
such a representation we analyze the magnetic field symme- 
try of the dc current flowing through the adiabatically driven 
scatterer in Sec lIIII We show that in the two terminal case 
there is a dc current 7' od ' that is odd in magnetic field, linear 
in lu and dependent on the applied voltage. To calculate cor- 
rectly /' od ' it is necessary to find the residual Floquet matrix 
A. Using several simple examples we outline the method for 
calculating A in Sec lIVI We conclude in SecfVl 



II. GENERAL APPROACH 

We use the scattering matrix approach-*— *— which views 
the mesoscopic sample as a scatterer which causes transmis- 
sion and reflection of incident carriers. The scatterer is as- 
sumed to be coupled to N r reservoirs via single channel bal- 
listic leads which we will number by the Greek letters a, f3, 
etc. 

We assume that in the stationary case electrons coming 
from the reservoirs and interacting with the scatterer are sub- 
ject only to elastic scattering. Such (single particle) scattering 
can be described with the help of the scattering matrix So- 
The index denotes the stationary scattering matrix. In gen- 
eral So is a function of the electron energy E. This matrix 
collects all the quantum mechanical amplitudes for electrons 
coming from some lead f3 to be scattered into the same or 
any other lead a. These amplitudes are normalized in such a 
way that their square define the corresponding particle fluxes 
(currents). If the electron velocities at a given energy are 
the same in all the leads we can use these amplitudes to re- 
late the incident and out-going wave functions. For instance, 

let ^0^0 (E> *) = e'^Vo!?^). be the amplitude of a wave 
function describing electrons with energy E incident in lead 
f3. Then the amplitude of the wave function of particles out- 
going in lead a, (E,t) = e _i # ^^(E), is defined as 
follows: 

4T ) (E) = T,s o , a0 (e)4:;\e). (i) 

0=1 

Current conservation implies that the scattering matrix is 
a unitary matrixi— 

SlSo = SoSl = I, (2) 

where / is a unit matrix. In fact, the knowledge of the matrix 
So(E) is equivalent to the knowledge of the solution for the 
stationary Schrodinger equation. 



For the dynamical problem with time-dependent scattering, 
scattering is characterized by the integral scattering operator 
which depends on two times—. One time-argument relates 
to the incoming states and the second time-argument to the 
outgoing states. In this paper we are dealing with a partic- 
ular non stationary case, namely with a periodically driven 
scattering problem. We assume that the scattering potential 
(hence the scattering properties of a sample) is varied in time 
periodically with period T = 27r/c£. Then, according to the 
Floquet theorem (see, e.g., Refs. |48h52T) . the solution for the 
time-dependent Schrodinger equation can be represented in a 
relatively simple form 

oo 

<f{E,t)=e- l i t e- inwt i>{E n ). (3) 

n= — oo 

Here E is the Floquet energy; ip(E„) is a general solution 
of the stationary Schrodinger equation corresponding to the 
energy E n = E + nhu. 

Scattering on such an oscillatory scatterer can be described 
via the Floquet scattering matrix. In this work we are con- 
cerned with the low-frequency properties of this dynamic 
problem and the relevant Floquet matrix Sf describes the 
transitions between the propagating states onlyp— The ele- 
ments SF, a 0(E n , E) of this matrix are the quantum mechan- 
ical amplitudes (normalized for current) for an electron with 
energy E to enter the scatterer through lead j3 and to leave 
the scatterer with energy E n = E + nhui through lead a. 

In particular, if the reservoirs are stationary then the in- 
coming wave function is 9 p(E, t) and the wave function for 
particles outgoing to lead a is of the form Eq.JHJ with 

y4 out) (£n) = E E \/ ^S F , a 0(E n ,E m )^(E m ) (4) 



Here k n = \/2m e E rl /Ti with m e being the electron mass. 
Physically Eq.@ means that an electron interacting with an 
oscillating scatterer can gain or lose one or several energy 
quanta nhu, n — 0,±1,±2, and thus an electron can 
change its energy by a discrete amount nhio. 

Current conservation implies again that also the matrix 
Sf is unitary. For the Floquet scattering matrix the analog 
of Eq.^ reads as follows: 

/J S F ,a0(En, E)SF,a~/{E n , E m ) = <5 m Q<5/3 7 , (5a) 



Sp i0l p{E, E n )SF,i0(E m , E n ) — S m o8 al . (5b) 

n 

Here the summation over n goes only over those n which 
correspond to a positive E n = E + nhui. In the low frequency 
limit we have Tiu> <C E, and thus n extends from — oo to +00. 

To find the Floquet scattering matrix one needs to solve a 
fully time-dependent Schrodinger equation. Compared to the 
stationary problem, this is a more difficult and, generally, it 
can be done only numerically. On the other hand the repre- 
sentation Eq.© seems effectively to reduce the periodically 
driven case to the stationary one. Therefore it is attractive 
to try to relate the Floquet scattering matrix Sf to the sta- 
tionary scattering matrix So. 
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A. Adiabatic approximation 

Let the stationary scattering matrix So(E, {p}) depend on 
a set of parameters pi g {p}, i — 1,2, ... , N p (e.g., the sam- 
ple's shape, the strength of coupling to leads, the magnetic 
field, etc.) . Varying these parameters one can change the scat- 
tering properties of a sample. We take these parameters to 
be periodic functions in time: p-i(t) = pi(t + T),Vi. Then the 
matrix So becomes time-dependent: So(E,t) = So(E, {p(t)}). 
In general the matrix So{t) does not describe the scattering 
of electrons by a time-dependent scatterer: only the Floquet 
scattering matrix Sf does. Nevertheless in the low frequency 
limit, w — *• 0, there exists a connection between these two 
matrices. This connection becomes more evident if one rep- 
resents the Floquet scattering matrix elements as a series in 
powers of u). 



1. zeroth order approximation 

To zero-th order in the driving frequency the elements 
of the Floquet scattering matrix Sf{E u ,E) can be approx- 
imated by the Fourier coefficients So,n of the stationary scat- 



tering matrix So as followsi— 

S F (E n ,E) = SoAE) + 0(u). 



(6a) 



S F (E,E n )=S ,- n (E)+0(oj). (6b) 

Here 0(uj) denotes the rest which is at least first order in 
frequency uj and which is neglected in the zero-th order adia- 
batic approximation. The Fourier transformation used reads 
as follows 



So(E,t)= e- m ^S , n (E), 



I 



dt 

r 



*S (£,t). 



(7a) 



(7b) 



Before proceeding we check that this approximation is con- 
sistent with the current conservation condition. Substituting 
Eq.© into Eq.JSJ and performing the inverse Fourier trans- 
formation we arrive at Eq.JS). 

Equation © corresponds to the frozen scattering matrix 
approximation. Within this approximation the stationary 
scattering matrix (with parameters dependent on time) com- 
pletely characterizes the time-dependent scattering. This ap- 
proximation is exact if the scattering matrix So is independent 
of the electron energy E within the relevant energy intervali— 



2. first order approximation 

To first order in the pump-frequency uj we can represent the 
Floquet matrix with the help of the frozen scattering matrix, 
its energy derivatives and a matrix A. In general the matrix 
A can not be expressed in terms of the stationary scattering 
matrix So and it has to be calculated (like So itself) in each 
particular case. The advantage of the representation which 
we introduce is that the matrix A has a much smaller number 



of elements than the Floquet scattering matrix. The matrix 
A depends on only one energy, E, and therefore it has N r x N r 
elements like the stationary scattering matrix So. In contrast, 
the Floquet scattering matrix Sf depends on two energies, E 
and E n — E+nhuj, and therefore has ~ {2n max + 1) X N r x ir- 
relevant elements. Here n max is the maximum number of en- 
ergy quanta Ttw absorbed/emitted by an electron interacting 
with the scatterer which we should take into account to cor- 
rectly describe the scattering process. For small amplitude 
driving we have n max ~ 1. In contrast, if the parameters 
vary with a large amplitude then n max 3> 1. We represent 
the Floquet matrix in the formi— 

S F (E n ,E) = S . n (E) + nf ™ dS °^ E) + hujA n (E) + 0(uj 2 ), 

Z UEj 



s F (E,E n ) = s ,-4E)+^ dSo -:" iE) +hujA„ n (E)+o(uj 2 ). 

z ohj 

(8b) 

Note the right hand side (RHS) of Eq. lSat is defined with re- 
spect to the incoming energy of carriers, while in Ea. ljSbfl the 
RHS is expressed in terms of the energy of outgoing particles. 
To first order in uj, the case of interest here, these two repre- 
sentations are fully consistent. Going from one representation 
to the other, one needs to take into account that the contri- 
bution from the first term on the RHS depends on the choice 
of the reference energy. The second and the third terms being 
themselves proportional to uj do not depend on this choice. 

In Eq.© we have introduced a new matrix A[E, i) with 
Fourier coefficients A n (E). The current conservation condi- 
tion, Eq.lJ^J, leads to the following equation for the matrix 



huj 



(sl{E,t)A{E,t) + A\E,t)S Q {E,tj) = 



V{Sl;S } = ih 



dSl 8S 



dt dE 



d§l dS 
dE dt 



(9a) 



(9b) 



Note the matrix "P{Sq;So} is traceless. Another but equiva- 
lent representation can be obtained from Eq. (l9at multiplying 
both sides from the left by So and from the right by Sq, and 
by taking into account that because of the unitarity condition, 
Eq.@, we have S d[Sl] = -d[S ]5^. 

We remark that Eq.@ tells us that the expansion in pow- 
ers of uj is, in fact, an expansion in powers of huj/8E, where 
5E is the energy scale over which the scattering matrix So(E) 
changes significantly. Therefore, the frequency uj can be con- 
sidered as slow and the expansion Eq.JHJ can be relevant if 

Tiuj < SE. (10) 
Consequently, to characterize scattering with an accuracy of 
order uj one needs to determine the matrix A. Equation © 
defines only the anti commutator of two matrices, So and A, 
and it is insufficient to determine the matrix A. 

By analogy with Eq.lJSJ we can express the Floquet scat- 
tering matrix elements up to first order in driving frequency 
in terms of the Fourier coefficients of some effective matrix. 
We introduce two matrices Si n and Sout defined with respect 
to incoming and outgoing energies, respectively: 



S in (E,t) = So (E,t) + 



m a 2 So 

2 dtdE 



+ hujA(E,t). 



(Ha) 
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S out {E,t) =S (E,t) 



ih d 2 S 



+ huiA(E,t). 



(lib) 



2 dtdE 

Performing the Fourier transformation of Eqs. llH and com- 
paring the result with Eqs.© we find: 



S F (E n ,E) = S in ,„(E) + 0(u 2 ). 



S F {E, En) = S out ,-n{E) + 0{uj 2 ). 



(12a) 



(12b) 



We emphasize that the matrices Si„(t) and S ou t{t) are not 
scattering matrices because they are not unitary: Their 
Fourier coefficients just define the corresponding matrix ele- 
ments of the Floquet scattering matrix according to Ea, l|12^ , 
Nevertheless these matrices conserve the current "on aver- 
age", i.e. after integrating over the time period T: 



I 



^Sl(E,t)S in (E,t)=i+0(u 2 ) 



(13a) 



^Sl ut (E,t)Sout(E,t) = I + 0(cu 2 ). (13b) 



Now we use Eq.© to analyze the general properties of the 
matrix A which are due to the micro reversibility of the 
Schrodinger equation with a periodically oscillating potential. 



B. Micro-reversibility and magnetic field 
symmetry of the Floquet scattering matrix 

We start with the stationary case when the single parti- 
cle Hamiltonian (and correspondingly the scattering matrix) 
is independent of time and recall some properties of the sta- 
tionary scattering matrix^-*— 

The micro-reversibility of the equation of motion (i.e., the 
Schrodinger equation) puts some constraints onto the scat- 
tering matrix. To make the notation more convenient let us 
arrange the incoming/outgoing wave functions at all the leads 
into the vector column 



/ in 



\ 



(14) 



\4>N r ) 

Then Eq.JTJ can be written in the compact form: 



(out) 



= S ijj 



(in) 



(15) 



The micro-reversibility condition (i.e., the invariance with re- 
spect to the time inversion) for the spin less case under consid- 
eration leaves the solution of the scattering problem invariant 
under the simultaneous inversion of the direction of move- 
ment, the inversion of a possibly present magnetic field H, 
and the replacement $ — ► \E r *. Therefore, the evolution of the 
two wave functions, namely ^{E,H,t) and ^* (E, —H, —t), 
is exactly the same and is described by the same scattering 
matrix So. Taking into account that the inversion of the di- 
rection of movement turns the outgoing waves to incoming 



ones and vice versa we can write the following equations for 
the starting solution and its transform: 



(out) 



(E,H) = S (E,H)^"\E,H), 



4> im) (E,-H)) = S (E,H)(ip {out >(E 1 -H) 



Uout), 



(16a) 



(16b) 



From the unitarity condition Eq.J^J it follows that S 1 = Sq. 
Therefore we can rewrite Eq. (ll6at as follows: ^ m \E,H) = 
Sl(E,H)4> {out) (E,H). Comparing the last with Eq.^TSEJ we 
arrive at the required condition!— 



(17) 



where the upper index "T" denotes transposition. 

Next we consider a periodically driven scattering problem. 
As we saw micro-reversibility requires the scattering matrix to 
be symmetric with respect to the interchange of incoming and 
outgoing channels. For the Floquet scattering matrix these 
channels are characterized by both the lead index and the 
number n showing how many energy quanta %u> an electron 
absorbs/emits during the scattering process. In addition, to 
get the required symmetry condition, we have to take into 
account that the parameters pi of the Hamiltonian depend on 
time. We suppose they change periodically in time with the 
same frequency u, and with possible relative phase shifts Lpc 

Pi(t) = Pi,o +Pi,i cos(ojt + ipi). (18) 

In such a case time reversal implies the inversion of the sign 
of all the phase shifts tpi. Therefore, the Floquet scatter- 
ing matrix elements are subject to the following fundamental 
symmetry: 

SF, a f>(E,E n ;H,<p) = S Fi( , a (E n ,E;-H,-<p), (19a) 

or in a matrix form 

S F (E;-H,- t p) = S^{E;H,<p). (19b) 

Here E is the Floquet energy [see, Eq.(|5J]; tp denotes the set 
of all the tpi . 

Next we derive the symmetry conditions for the matrix A 
entering Eq.JHJ. Our definition of the phases ifii [see, Eg. 1181 1 
implies that the frozen scattering matrix So(E, i) [i.e., the sta- 
tionary scattering matrix with parameters dependent on time 
So(E,t) — So(E,pi(t)) ] possesses the following symmetry 

S (E,-t;H,- l p)=S (E,t;H, l p). (20) 

Then equation @ gives us: 

A(E,-t;H,-<p) = -A(E,t;H,<p). (21) 

Correspondingly, for the Fourier coefficients, we have the fol- 
lowing: 

S , n (E;H,-(p) = S ,- n (E;H,<p), (22a) 

An(E;H,-tp) = -A-n(E;H,<p). (22b) 

Substituting the equations given above into the adiabatic 
expansion, Eqs.©, and taking into account the micro- 
reversibility condition, Eq. (ll9ll . we find the required symme- 
try condition for the matrix A(t): 



A{-H) = -A T (H). 



(23) 
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In particular, in the absence of magnetic fields, H = 0, the 
diagonal elements of A vanish. That was previously shown in 
Ref. IHS Alternatively equation I23H can be obtained directly 
from Eq.@ exploiting the symmetry condition Ea. l|17|l and 
the unitarity of the frozen scattering matrix So(E,t). 

The symmetry properties of the residual Floquet matrix A 
are completely different from that of the stationary scatter- 
ing matrix So- The residual Floquet matrix A reflects directly 
the most important differences between an adiabatic scatter- 
ing process at a periodically evolving scatterer and a strictly 
stationary scattering process. 



III. MAGNETIC FIELD SYMMETRY OF THE 
DC CURRENT FLOWING THROUGH THE 
SLOWLY DRIVEN SCATTERER 

Now we use the results of the previous section to ana- 
lyze the dc current through the mesoscopic sample with pe- 
riodically varying parameters. We will consider two mech- 
anisms which can give rise to such a current. The first 
mechanism is a quantum pump effect consisting in rectifying 
of time-dependent currents generated by the non stationary 
scatterer?— Second we permit a constant in time difference 
of electrochemical potentials/temperatures between the dif- 
ferent reservoirs. The last is important, because the widely 
investigated situation with reservoirs being at the same elec- 
trochemical potential actually hides some physics underlying 
the quantum pump effect. 

The dc current I a flowing from the scatterer to the reservoir 
in the lead a can be calculated as follows!— 

I a = l JdE \S F MEn,E)\ 2 f 0i p(E) ~ / , Q (£) j . 

(24) 

Here fo, a is the electron distribution function for the reservoir 
a. We assume that the reservoirs are in a stationary equilib- 
rium state with possibly different electrochemical potentials 
(j, a and temperatures T a . Then fo. a is the Fermi distribution 
function 

foAE) = ^=^, (25) 

1 + e h B T a 

with fcfl being the Boltzmann constant. Substituting the adi- 
abatic expansion Eq.JSJ into Eg. 12411 and performing the in- 
verse Fourier transformation we find the current up to linear 
in uj terms as follows!— 



00 t , 
I a = fdEf$ E\foAE) 

n n a l - 



dI af3 (E,t) 
dE 



(26) 



+ j;\S , al 3(E,t)\ 2 [f O:0 (E) - /„,„ (E)] } , 



where dI a p/dE is a spectral current driven by the non sta- 
tionary scatterer from lead (5 into lead a: 



dIaf3 = ~ ( 2huRe[Sl aP A a p] + ^V{S , al3 ; S^} ) . (27) 



dE 



Here 7?e[X] is the real part of X; the function "P{X;Y} is 
defined in Eq. l9bll . The spectral currents dl a ^jdE are subject 
to the following conservation lawi— 



E 



dI a0 (E,t) 
dE 



0. 



(28) 



Using Eq. I28H and the unitarity of the frozen scattering ma- 
trix, J2 a \ S o,ap\ 2 = J2p \ s o, a p\ 2 = 1, one can easily check 
that the current I a is conserved: ^ Q I a — 0. Further, using 
the symmetry conditions, Eqs. 11171 and I23II . and rearrang- 
ing the terms in Eg. 1261 we divide the current into the even, 
li ev) (H) = li m \-H), and odd, li° d) (H) = -l£ d \-H), in 
magnetic field parts: 

00 T ( 

/W(j?) = f/ffi/fE I/o^-W 
1 



+ hu}Re[So, a0 A a/3 - So :0a A0 a ]) (29a 



+ L/0./3 + /o, Q J 1 - e -j 



00 T ( 

ii od) (H) = %fdEf§j:\[f , g -f . a ] 

( 
./IfiwP-ISo^ , P{So,^;So, Q/3 }-y{So,gc,;S ' gc ,} 

1 2 1 4 

+[fo,p + fo, a ]huRe[S^ al3 A a p + So^Ap 

(29b) 

To show that both these currents are separately conserved, 
i.e., that J^Jq ' = and E Ia d ^ = 0, one can use the 

a a 

relations— 



AUuj J2 Re[Sl a0 A al3 ] = V{Sl So} 



BP, 



(30a) 



4hto £ Re[So, al3 A af3 ] = P{S ; Sl} aa . (30b) 

0=1 

which follow from Eo. l|9a|l . 

In a general multi-terminal situation, i.e., if not all the 
reservoirs are at the same potential (temperature), the main 
contributions to both the even li e "' and the odd li 0<i ' currents 
are proportional to the conductances ^-|iSo,a/3| 2 averaged over 
time. The non-stationarity results only in small corrections. 
However in the two terminal case the odd in magnetic field 
dc current /( od ' has no contribution coming from the conduc- 
tances. The current /' od ' is linear in uj and it is entirely due 
to the non-adiabaticity of the pump scattering processes. 



A. Two terminal many channel scatterer 

To show this let us consider the scatterer connected to only 
two reservoirs via, possibly many channel, ballistic leads. We 
will mark the quantities related to the left and to the right 
reservoirs via the lower indices " L" and " R" , respectively. Let 
the left lead have Nl channels, and the right lead have Nr 
channels: Nl + Nr — N r . We define the currents flowing to 
the left II and to the right Ir = —II, and the distribution 
functions for the left /o,l and for the right fo,R reservoirs as 
follows: 



E 1° 

Q = l 



fo,a = fo,L, 1 < a < N L , 



Ir= E 7 «, fo, a = fo,R, N L + 1 < a < N r 

a=N L +l 
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By analogy we redefine the quantities dependent on two in- 
dices. For example, the reflection to the left Rll and the 
spectral current (Url /dE driven from the left to the right are 
defined as follows: 

N L N L 

Rll = E E \So,np\ , 

Q = l = 1 



dE 



Nl 



2-/ dE • 



c=n l +i g=i 



Note that the two terminal transmission is symmetric in reser- 
voirs indices, Tlr = Trl, and it is even in magnetic field. 
That can be easily seen from their definition, similar to the 
one given above for Rll, and from the unitarity of the scat- 
tering matrix So, Eq.@. In addition from Eo. 1281 we get: 
dI Li jdE + dI R ^/dE = 0, for £ = L, R. Using the identity^ 



is odd in both V and AT. Both contributions, j( od -™) and 
j(od,od) ^ nave -j-jjg same origin: They are rectified ac currents 
with spectral density dI a p/dE, Eg. 12711 . pushed by the pump 
from one reservoir to another. The part /( od ' e ") emphasizes 
the contribution arising if there are incoming electrons from 
both leads, a and /3. While the part j r ( od > od ) ; s entirely due to 
an asymmetry in electron flows incident from the leads. This 
asymmetry, due to the difference between the reservoir's dis- 
tribution functions /o, Q and fo,g, vanishes in the absence of an 
applied voltage, V = 0, and in the absence of a temperature 
difference, AT = 0. 

For further reference we now give the equations 1311 for the 
particular case of a scatterer connected to one-channel leads. 



1. Two terminal single channel scatterer 



dI L L dI L R 



dE 



dE 



d! L 
dE 



2tt 



dSodSl _ dS dSl 
dt dE dE dt 



performing necessary summations in Eas. 1291 . and integrat- 
ing by parts over time and over energy, we get: 



For single channel leads, Nl = Nr = 1, the stationary 
scattering matrix So is a unitary 2x2 matrix: 



So 



^Re~ ie i 
is/Te i4 



VTe- i4 
VRe w 



(33) 



T (ev) 



JdEj&\[f , R -f , L ] 



x(tlr+j: e 

a=10=N L +l 



dt 

h J J T 



V y{gQ,og;g;, o g}+y{5 0- g c ,;S * /3Q } 



+ I 



a 

" dE 



[fo,R + /o,l]) f E (#St + St%) 

a=l v ' <■ 



(31a) 



r(od) 



r(od,ev) | r(od,od) 
1 L ~^~ 1 L ' 



j-(od,ev) 



r(od.od) 



OO 1 r 

§JdEj^{(-^[f ,n + fo, L ]) 

V ( SSo at St SS \ \ 

2^ l~ST°o s q — ) h 

v = 1 \ / OLOL s 

oo T , 

JdEjf{[f ,R-fo, L ] 



(31b) 



2vr 



X E E Re[So, a p A *0 + SZ t p a A0 a ] >■ 

a=10=N L +l J 

For low driving frequencies, lu — > 0, we see that in the two 
terminal case the part of the dc-current that is odd in mag- 
netic field, I {od) (H) = -^(-H), is linear in u irrespec- 
tive of whether the reservoirs are at the same conditions 
(fo.L = fo,n) or not (f , L / fo,n)- 

Let us introduce the voltage V and the temperature differ- 
ence AT applied to the system, both constant in time: 



/J-r = Mo + ^f, 



T R = To + 



eV 

U>L = UQ —, 



T L = To 



AT 
2 ' 



(32) 



and analyze the current /( od ) in more detail. According to 
Ea. )|31b0 this current consists of two parts, I 



r(od,od) 



(od) 
L 



*(od,ev) 



The first one 
I ( ° d - ev \V, AT) 



r(od,ev) 



(-V,-AT), 



is even in both V and AT and it survives even at /o,l = 
fo,R- This contribution is due to conventional quantum pump 
effect. 8 In contrast, the second part 

I ( L od ' od \V,AT) = -I { ° d - od \-V,~AT), 



Here R and T are the reflection and the transmission prob- 
ability, respectively (R + T = 1). The phase 6 characterizes 
the asymmetry between the reflection to the left and to the 
right. The phase 7 relates to the charge on the scatterer. The 
phase (j> characterizes the asymmetry between the transmis- 
sion through the scatterer from the left to the right and back 
and it relates to the magnetic flux on the scatterer. 

We assume that all these quantities are functions of the 
electron energy E, the magnetic field H, and the external 
parameters Pi(t) varying with frequency uj. From Eq. 1171 it 
follows that R, T, 7, and 6 are even functions of the magnetic 
field H , while <f> is an odd function of H. 

Using the scattering matrix, Eg. 1331 . we rewrite the dc 
current II = — Ir, Eg. 1311 . as follows: 



j{ev ,ev) 



+ 1 



(ev,od) 



+ 1 



(od,ev) 



+ 1 



(od,od) 



(34a) 



00 T 

j-{ev,ev) & f t r- , f dt ( d 



■ -to I ""J T\-8E^* + f °^ ) R W (34b) 



06 



00 T 

i*-~-5/«/4W-m(t-!£[t£] 



(34c) 



00 T 







l( od,od) = ^_J dE J ^[f 0iR -f , L ]R e [S* , LR ALR+Sl RL ARL] 


(34e) 

Here the first upper index, ev/od, relates to the magnetic 
field symmetry of the current, while the second upper index 
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relates to the symmetry with respect to the applied voltage 
(temperature) difference. 

The currents j^ v ' ev > ant [ j(° d ^ v ) are conv entional pumped 
currents (with reservoirs being at the same conditions) . They 
depend on the asymmetry of the stationary scattering matrix: 
The phases 9 and <f> describe the asymmetry in the reflection 
from and in the transmission through the scatterer, respec- 
tively. 

The remaining two contributions, i ( £ v > od ) anc j j(od,od) ^ are 
present if the electron flows incoming from the reservoirs are 



(ev ,od) 



exists 



different. The even in magnetic field current I 
already in the stationary case. The variation of the scattering 
parameters results in averaging over the time period and in 
the correction to the frozen conductance. 26 In contrast the 
odd in magnetic field current /^ od ' oci > exists only in the non- 
stationary regime. 



{ev ,ev) j(od,ev) j(od,od) 



The contributions I 
(proportional to dy/dt) of / 



(ev,od) 



I L ' , and the part 
all are due to the quantum 



rectification of ac currents, Eg. 12711 . generated by the oscillat- 
ing scatterer. This mechanism does work (i.e., a dc current 
exists) if the time reversal invariance is broken in the system 
by the varying parameters pi and hence the integral over the 
time period does not vanish. 

The last statement is not evident for the current j(° d ' od \ 
To make it clear we note the following: Generally the non- 
adiabatic corrections to the frozen solution of the Schrodinger 
equation (and to the frozen scattering matrix) are propor- 
tional to the time derivative. The matrix A being a part 
of these corrections should be proportional to d/dt as well. 
Therefore, the conditions for ' od ' ; Eq. l34eH . to be non- 
vanishing, (which requires that the integral over the time 
period is nonzero), are generally the same as that for, e.g. 



Eq.JM 

However, we emphasize that the necessary conditions to get 
a pumped current include both the time reversal symmetry 
breaking and the presence of a spatial asymmetry. Strictly 
speaking these conditions are not identical for all the parts, 
Eqs. H34bll - I34el . In particular, the applied voltage V makes 
the whole system (the sample plus reservoirs) to be spatially 
non-symmetric, in the sense that the direction from the left 
to the right is not identical to the opposite direction. There- 



(od,od) 



is less sensitive to the spatial 

(od,ev) . 



fore, e.g., the current I 

asymmetry of a scatterer than, e.g., the current I 

In the next section we calculate the residual Floquet ma- 
trix A for several simple examples and illustrate the validity of 
the general statements made here. In particular we consider a 
one-dimensional loop with two leads and with enclosed mag- 
netic flux. This example shows that the current j(° d '° d ^ exists 



(at nonzero magnetic flux) already if only the time reversal 
invariance is broken (i.e., if two parameters oscillate with a 
phase lag Aip 7^ 0). In contrast, a non-zero current j(° d ' ev > 
requires (in addition to Aip 7^ 0) an asymmetry in coupling 
to the leads. If such a coupling is symmetric then the current 
j(od,ev) j g identically zero no matter how much the magnetic 
flux through the ring is and whether Aip is zero or not. 



IV. THE FLOQUET SCATTERING MATRIX 
FOR LOW DRIVING FREQUENCIES: SIMPLE 
EXAMPLES 

In this section we illustrate how one can calculate the linear 
in lu corrections to the frozen scattering matrix in the same 
fashion as the stationary scattering matrix So- 

According to Eg . 11121 1 . at u — > the Floquet scattering 
matrix elements are the Fourier coefficients of some matrices 
Sin/Sout, Ea. dlH which depend on the stationary scattering 
matrix So and on the matrix A. The matrix Si n / S ou t does 
not possess a definite symmetry, i.e., with respect to a time 
and/or a magnetic field direction reversal. While the sta- 
tionary scattering matrix So and the matrix A do. The last 
circumstance is the motivation why we expressed the current 
I a in terms of So and A instead of Si n /S ut- On the other 
hand for the calculation of the Floquet scattering matrix el- 
ements it is more convenient to work in terms of the matrix 

Sin Or Sout- 



A. Single 5-function barrier 

In the first example we consider the Floquet scattering ma- 
trix in the limit ui — *• of an oscillating point-like scatterer 
coupled to two reservoirs via one-channel leads. As we will 
show for such a scatterer 



,4 = 0, 



(35) 



and low frequency scattering up to linear in w terms is entirely 
described by the frozen scattering matrix So(t), see Eqs. ll2H 
and (TTJ at A = 0. 

To find Sf we have to solve the Schrodinger equation with 
the potential V(x, t) being the delta function S(x) multiplied 
by the amplitude oscillating in time: 



V(x, t) = 8{x) (Vo + 2Vi cos(wi + ipfj . 



(36) 



According to the Floquet theorem the solution of the above 
equation has the form of Eq.^. Away from the point x — 
the functions ip(E n ) are the plain waves: 



ip(E n 



(37) 



The coefficients a n , b n are determined from the boundary 
condition at x = 0: 



*(x = +0) = 9(x = -0), 



a* I 

dx lx=+0 



8* 
dx I 



2g?-V(t)*(x = 0). 



(38) 



First, to find Sf,ll and Sf.rl we consider the plain wave 
of a unit amplitude with energy E coming from the left (we 
directed the x-axis from the left to the right): 



* (m) (£,t) 



(39) 



Here E = fi k /(2m e ). Then the coefficients a„ u ' and b 
for an outgoing wave 



, (out) 



(out) 



(x)a { ° wt) e iknX + 1 



-x)b { ° ut) e- iknX 



(40) 
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(here 8(x) is the Heaviside step function: 6(x) = 1 at x > 
and 8(x) = at x < 0) define the Floquet scattering matrix 
elements as follows: 



SF,RL{E n , E) = \ 

Sf,ll(E u , E) = J -ff-b 



(out) 



(out) 



(41) 



Substituting the whole wave function * = * l * nJ + ■^ {out> into 
the boundary condition Ea. l|38p we get the following relations 
between the different al? ut ^ and bn Ut ^\ 



(k n + i/t )aL° ut) = k n 8 n0 - i(Kia { °^ + «_ 1 afc 1 ^), 



, (out) 
fJ-n 



(out) 



(42) 



SnO- 



Here we have introduced the following parameters: 



Alternatively, one can use Eq.© to show that the ma- 
trix A vanishes for the oscillating <5-function potential. It 
is because the commutator V{§l;So} is identically zero for 
the scattering matrix Eg, 1471 . We can conclude that a point 
scatterer can not generate a quantum pump effect since it can 
not rectify ac currents (the spectral density dI a p/dE, Eq. l27L 
vanishes). An oscillating scatterer does of course generate ac- 
currents, but these currents are total time derivatives of the 
charge near the barriej— and thus can not contribute to a 
dc- current. 

Note, that the deviation of the effective scattering matrix 
Si„(E,t), Eg. 1461 from the frozen scattering matrix So(E,t), 
Eg. 1471 . is as small as, at least, Tiw/E. For the opaque barrier 
the deviation is even smaller due to the factor k/n <C 1. For 
the small oscillating amplitude case the deviation is addition- 
ally damped by the factor kx/kq <C 1. 



-Vo, K±i 



-V x e 



(43) 



We solve Eo, l|42^ in the adiabatic limit u> — > of interest 
here. In this limit we can safely expand the wave vector k„ 
as follows: 

k n = k+ — + O(co 2 ), (44) 



where v = hk/m e is an electron velocity. In addition we use 
the adiabatic expansion Eq. H12at and express Sf{E u ^ E) in 
terms of the Fourier coefficients of the matrix Si n (E). Sub- 
stituting Eqs.(0^, pijl and l|12a[l into Ea. lH^ . and ignoring 
all the terms of order u 2 and higher we can write: 

(k + iKo)Sin,RL,n + (j — *ffe) ^So,RL,n = 
k5n0 — i(KlSin y RL,n-l + K- 1 Sin,RL,n+ 1 ) 
+ «2^fc [tl(n — l)So,RL,n-l + K-l(n + l)S(),RL,n+l], 



B. Scatterer composed of two point-like barriers 

In this subsection we consider an example of a spatially 
"extended" scatterer which consists of two point-like scatter- 
ers placed at x — and x = L, respectively. This system is 
coupled to two reservoirs via single channel leads. 

The scattering properties of point scatterers are assumed to 
be oscillating in time with the same frequency ui. Scattering 
at the left and on the right barriers is described via the Flo- 
quet scattering matrices Sp and Sp, respectively. Scattering 
on the whole system is described via the Floquet scattering 
matrix Sf. 

By analogy with the previous example we consider scatter- 
ing of a unit wave coming from the left, Eg. 13911 . The whole 
wave function is of a Floquet function type Eq.JjU with 



Si„ t LL,n(E) — Sin,RL,n{E) — 8nO\J ^fr ■ 

Performing the inverse Fourier transformation we find the 
equation for the time-dependent matrix elements of the ma- 
trix Si n (E, t): 



Si n ,RL(E, t) 



i k-iK(t) dS , RL (E,t) 
2vk k+in(t) dt ■ 



k + in(t) 

Sin,Lh{E,t) = Sin,RL(E,t) — 1 



(45) 



Here K,{t) = m e V{t)/h 2 . We solve these equations perturba- 
tively in the small parameter proportional to d/dt ~ uo — > 0. 
To find the matrix elements Si„ t RR and Si n ,LR one can either 
exploit the symmetry condition or solve the same problem but 
with the unit wave incoming from the right: \]A m ' (_E, t) — 
e~ 1 ~K t e~' kx . Up to terms linear in uj, the solution of both 
problems reads: 



S in (E,t) = S (E,t) + 



ih d 2 S (E,t) 
2 dtdE ■ 



(46) 



Here we used dk/dE — l/(hv). The stationary matrix is well 
known: 

Comparing equations 146H and lllat we arrive at the an- 
nounced result, Eq, I35H . Thus, to describe the low frequency 
scattering on point-like scatterer it is enough to know only 
the frozen scattering matrix. 



n,oe 



i>(E n ) = I a n e lk ^ + b n e~^ x , < x < L, . 

x > L, 

(48) 

To find the unknown coefficients we use the boundary condi- 
tions which we formulate in terms of scattering matrices Sp 



+ \/i : S F , LL (E n ,E)e- 

-ik n x 



Sf,ll(E?zi E) = S FjLL (E n ,E) + J2 Sf,lr(Eu, E n 

m 

T fl n = Sp iRL (E n ,E) + J2 S F:RR (E n ,E m )J l^fc. 



-b n e 



= J2 Sp )IjL (E n ,E„ 




Sf,Rl(Eti, E) — Sp.fLL^En, Em) \ p-dmC 



(49) 

To simplify this system of equations we use the adiabatic ap- 
proximation Eq. 1121 for the Floquet scattering matrices. For 
this approximation to be valid, the energy quantum hui should 
be small compared with the relevant energy scale for the prob- 
lem, [see, Eq.©]. 

In the case under consideration, there are several energy 
scales. The first one is determined by the energy E of an 
incoming electron. This scale relates to the deviation of the 
effective scattering matrices Sf^ and S^ for point-like scat- 
terers from the corresponding frozen ones. This deviation is 
of the order of Tiuj/E. Another energy scale SE relates to 
the spatial size of the system L and arises from the quantum 
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mechanical interference in the region between the scatterers 
at < x < L. In our case, Eo. l|49[l . the interference effect 



is described via the factors e 1 ' 
follows: 



which we will expand as 



(50) 



Here u>l = v/L defines the distance AE ~ Tilol between the 
quantum levels if the system is decoupled from the reservoirs. 
The second term in the brackets on the RHS of Eq. )fi 
is due to an interplay of a quantum-mechanical interference 
with a quantized energy exchange between the scatterer and 
an electron traversing it. 

The system can be treated as spatially "extended" if L 2> 
Xe, where Xe = h/ \f2m e E is the de Broglie wave length for 
an electron with energy E. In such a case the non-adiabatic 
corrections to the frozen scattering matrix are at least of or- 
der Tiuj/AE. Note if the energy E is close to the energy 
of a transmission resonance then the corrections will be of 
order Tuj/T, where F is the width of the transmission res- 
onance. In contrast, if L <C Xe then the scatterer can be 
viewed as point-like and the non-adiabatic corrections will 
be as small as hco/E (see, Sec II V Al . Therefore, assuming 
I > we can safely ignore the corrections of order hio/E 
and concentrate on the larger corrections of order hu/8E with 
SE = min{A_E,r}. Since we ignore the terms of order hu>/E 
we can replace the Floquet scattering matrices for point-like 
scatterers by the corresponding frozen scattering matrices, 
S R/L (E n ,E m ) = §o,L-m + 0(huj/E) [see, Eq.©]. To avoid 
a possible misunderstanding we do not write the energy E 
as an argument of S R ^ L emphasizing that these matrices can 
be treated as energy independent on the scale of order SE. 
Nevertheless they can depend on energy over a much larger 
scale, say, of order E. On the other hand, since we keep the 
terms of order hu>/8E we use the adiabatic approximation 
S F {E n ,E) = S in ,„(E) + 0(u 2 ) [see, Eq.(T2)] for the Floquet 
scattering matrix of the whole structure. 

Using these approximations and substituting Eq. l5UI into 
the system of equations 1491 and performing the inverse 
Fourier transformation we arrive at the following time- 
dependent equations valid up to first order in d/dt: 

S in , LL {E,t) = S^ LL (t) + S^ LR (t)b(t), 
(*)• 



o(t) = SoV 



st RR {t) bit), 



S ln ,R L (E,t) = S R RL (t)e' kL (a(t) 



1 da(t) 
wt. dt 1 ' 



1 da(t) 
TJ7 dt I ' 



(51) 

Here we introduced the functions a(t) and b(t) defined as fol- 
lows (x = a,b): 



(52) 



We consider the terms da/dt and db/dt as small perturbations 
and solve the system of equations (1511 up to linear order in 
these corrections terms. 

Note, that without the terms da/dt and db/dt the system 
of equations Eg. 15111 is exactly the system of equations which 



defines the matrix elements of the frozen (stationary) scatter- 
ing matrix (with the evident replacement Sin — * So)- 

Analogously, to calculate Si n ,RR and Si„,LR we consider the 
same problem but with the unit wave coming from the right: 
* (in) (_E,t) = e ~ % ~ t e~ ik{x ~ L) . It is convenient to represent 
the results in the matrix form: 



Sin(E, t) = So 



ujl at L 



Here So is the frozen scattering matrix: 

S (E,t) = M + MlM^Mr. 



(53) 



(54) 



The matrices M are all expressed in terms of the scattering 
matrix elements for the left and right scatterers. They depend 
on energy through the factor e' kL and on time through the 
matrices Sq and S R : 



M 



Ml 



°0,LL 







OR 

°0.RR 



M 



So.LR 



M R 



qR %KU Q 

°0,RL e u 

Our aim is to calculate the matrix 

TlUlA = Sin — So - 



1 st 

-S?LLe i2kL 1 



Sq,rl 

u °0,LR e 



(55) 



ih d 2 Sp 
2 dtdE 



[sec, Eg. lllall l for the double scatterer structure under con- 
sideration. Using Eg. 1541 we obtain 



m a 2 S 



1 d 

~2uj l dt 



Ml M 



2 Mr 



2 dtdE 

Then using Eg. 1531 we obtain 

hioA(E,t) = gi_ {A [MiM" 1 ] M-HiR 



(56) 



The advantage of the above expression is its compactness. 
However to be able to draw some conclusion we need the 
matrix elements of A expressed directly in terms of the matrix 
elements of So and S R : 



TiujAlr = 



So,LL-S a LL _d_ 

uj l a at 



In 




TiluAl 



,a !(2-A)& 



In 



In 



(57) 



TiojArr = 



So.rr-S RR Q 

U L A dt 
L qR A2kL 



hi 



Here A = 1 - S^ RR S R LL e rAKlj . We see that for the case 
considered the matrix elements of A are proportional to time 
derivatives as it should be. 
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FIG. 1: A one-channel ring of length L — L\ + Li with 
enclosed magnetic flux $ and with two leads. The Greek 
letter a numbers the scattering channels at the left S L and 
right S R wave splitters. 




0.4 0.6 

A(p, 2% 



1.0 



We can now use Eo. l|57[l to investigate the symmetry prop- 
erties of the matrix A. We suppose that the matrices 5^ and 
So* are of the form given by Eg. (1331 . If there is no magnetic 
field, H = 0, then in Eq.^ the phases <j> L/R = 0. There- 
fore, the scattering matrices are symmetric in the lead indices, 
Stp R — Sp^ R . In such a case A aa = and Alb. = — A RL 
that is in agreement with Eg. 1231 . In addition, if the system 
is inversion symmetric, i.e., if the scatterers are the same, 
S L = S R , and if in addition 8 = and = 0, then all matrix 
elements are zero, A = 0. Note, for an oscillating scatterer 
the last two properties apply only if two scatterers oscillate 
in synchronism, but not if they oscillate with a phase lag. 



FIG. 2: The parts of the pumped currents which are odd in 



magnetic flux i^° d ' ev ' ) 



and / 



(od,od) 



Eos. ll34dl and llSIel . 



given as a function of the phase difference Atp = <p — tp 
between the oscillating coupling parameters e R {t) and £ L (t) 
for symmetric coupling: Both the left and the right wave 
splitters are coupled symmetrically to the branches of the ring. 
The currents are in units of eu/(2n). The parameters are: 
Li = tt; L 2 = 0.7tt; $ = 0.4$ ; Mo = 9; V = 0.5; T L = Tr = 



0.1; e£ 



e R = 0.5; ef 1 = ef = 0.4. We use the following 



units: 2m e = 1; h = 1: e = 1: kn = 1. 



C. Ring with enclosed magnetic flux 

Our third example is a one-channel ring with enclosed mag- 
netic flux $ coupled to two reservoirs, FigQ The lower and 
the upper branches of the ring have length L\ and L2, respec- 
tively. Following Ref.|54|we describe the coupling between the 
ring and the lead via a single parameter 3x3 scattering ma- 
trix So(e), with e = e L and e = e R for the left and the right 
coupling points, respectively. The numbering of scattering 
channels is shown in Fig0 We will use two different ma- 
trices, Sq and Sq , which couple the lead to the branches 
of the ring symmetrically and asymmetrically, respectively. 
They are: 



-(a + b) 




y/e a b 


(58a) 


\fl b a J 




a yfi b \ 




-{a + b) 


(58b) 


b x/e a J 





Here a = (VI - 2e - l)/2 and b = (VI - 2e + l)/2. The 
coupling parameter should be within the following interval: 
< e < 0.5. We suppose that the coupling parameters e L 
and t R oscillate in time with the same frequency u) but with 



the phase lag Ay> = (p R 



L/R 



r L/R , L/R 1 , 
|e + e i cos(cjt - 



L/R 



)]/2- 



We keep the parameters e L and e R independent of the elec- 
tron energy. This allows us to describe the time-dependent 



scattering at the three lead splitters within the frozen scat- 
tering matrix approximation, see Sec llI A T1 Also we assume 
that a voltage V can be applied between the reservoirs, here 
kept at the same temperature, see Eg. 1321 . 

We concentrate mainly on the part of the current I L ° d ^ 
which is odd in magnetic flux. This current consists itself of 
two parts, 7 



(od,ev) 
L 



and I 



(od,od) 



, which are defined by Eas. l34dl 
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FIG. 3: The currents I { L od ' ev) and I L od ' od) [in units of ew/(27r)] 
as a function of the enclosed magnetic flux $ [in units of the 
magnetic flux quantum $0 = h/e] for the case of symmetric 
coupling: Both the left and the right wave splitters are cou- 
pled symmetrically to the branches of a ring. The parameters 
are the same as in Fig|2] The phase difference is: Aip = 7r/2. 
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FIG. 4: The currents through the ring are given as a function 
of the Fermi energy ^to- The coupling is asymmetric: Both the 
left and the right wave splitters are coupled asymmetrically 
to the branches of a ring. The parameters are the same as 
in Fig|5| The phase difference is: Aip — tt/2. The currents 



j(od,ev) 



and I 



(od,od) 



even in magnetic flux j^ v '° d ^ 
of 1/(2tt). 



are given in units of eu)/(2n). The current 



is given in dimensionless units 



and 134el l. respectively. To find J^° d > od ^ jt is necessary to 
calculate the residual Floquet matrix A for the system under 
consideration, Fig0 These calculations are quite similar to 
those given in Sec liVBl and we do not give the details here, 
see Appendix. 

The current j^° d > ev ^ odd m magnetic flux and even in ap- 
plied voltage, is governed by the phase <j> which determines 
the asymmetry in the transmission phase through the ring, 
see Ea. l|38|l . Interestingly, in our model this asymmetry de- 
pends crucially on the type of coupling between the leads and 
the ring. If each of the leads is coupled symmetrically to the 
arms of the ring, i.e., if = S ( s) (e L ) and S£ = S { s) (e R ), 
then the phase <f> = for any magnetic flux. In such a case 



the current / 



(od,ev) . 



is identically zero and the full current odd 



in magnetic flux Ij° d ^ is odd in the applied voltage as well. 

In Fig|2]and Fig|2]we depict the current contribution which 
is odd in magnetic flux as a function of the phase lag and 
the magnetic flux, respectively, for symmetric coupling. Thus 
if the coupling between the ring and the leads is symmet- 
ric the current odd in magnetic flux appears only at V 7^ 
(and/or at AT 7^ 0), i.e., when the electron flows incident on 
the ring from the reservoirs are different. In contrast, if any 
lead (or both) is coupled asymmetrically to the ring, i.e., if 

(e fl ), then the phase <f> is 
does contribute 



Ma) 



So = S^ a) (e L ) and/or §g = S' ( 
not identically zero and the current jj? ' ev > 



to I 



(od) 



In Fig|l] we depict the currents / 



(od,ev) 



and I 



(od,ev) 



function of the Fermi energy fio for asymmetric coupling. For 
comparison we give the current J^ v ' od '^ Eo, l34cL which is 
proportional to the transmission probability through the ring. 



Both currents I 



(od,ev) 



and j(° d ' od > peak (by modulo) at en- 
ergies where the transmission (reflection) coefficient changes 
sharply. 



V. CONCLUSION 

In this work we analyze the scattering properties of a peri- 
odically driven mesoscopic scatterer. Traversing such a scat- 
terer an electron can gain or lose one or several energy quanta 
TvjO and thus can change its energy. Therefore, generally the 
scattering matrix of a periodically driven mesoscopic scat- 
terer depends on two energies, incoming and outgoing, We 
show that at low driving frequency u — > one can introduce 
effective matrices depending on only one energy, either incom- 
ing or outgoing [see, Ea. Hll|l ]. which approximates accurately 
the Floquet scattering matrix up to terms of order u> [see, 
Eg. 1121 1. We introduce two effective matrices, Sin and Sout, 
which are not unitary. Nevertheless each of them conserves 
the current after averaging over a driving cycle. 

The matrices Si n and S ou t are the sum of a frozen scat- 
tering matrix and a matrix which determines the linear in u> 
part. The last is responsible for the quantum pump effect 8 
and it consists of two contributions. The first one is the sec- 
ond derivative of the frozen scattering matrix So(t). The sec- 
ond contribution is defined by an in principle new matrix A. 
In particular, the matrix A has a symmetry with respect to 
magnetic field reversal, Eg. 1231 . that is opposite to that of the 
stationary (frozen) scattering matrix, Eq.QJJ. In contrast to 
the stationary scattering matrix the residual Floquet matrix 
reflects directly the chirality of the pumping process. 

Using the adiabatic representation Ea. l|12|l for the Floquet 
scattering matrix we examine the dc current flowing through 
the two terminal (many channels) mesoscopic sample under 
the simultaneous action of a slow parametric oscillation of the 
scatterer and simultaneously subject to an applied dc volt- 
age. We divide the current into parts with definite symmetry 
properties with respect to a magnetic field and/or a voltage 
inversion. 

As it is known in the stationary case the dc current through 
the coherent two terminal sample is an even function of a 
magnetic field. On the other hand the periodically driven 
scatterer shows an odd in magnetic field, linear in u) current, 
Eo.pTb). which is due to the quantum pump effect. The odd 
in applied voltage part of this current is proportional to the 
residual Floquet matrix A [see also, Ea. l|34e|l ]. 

We demonstrate that the calculation of the residual matrix 
A can be performed in close analogy with the calculation of 
the stationary scattering matrix So- We emphasize that the 
matrix A reflects the interplay of absorbing/emitting of en- 
ergy quanta hu with quantum mechanical interference inside 
the scatterer. For instance, for a point-like scatterer (without 
the space for interference inside) the matrix A is identically 
zero. 

Our work suggests that additional experiments which in- 
vestigate a driven mesoscopic conductor in a less symmetric 
setup, i.e., with reservoirs having different electrochemical po- 
tentials or temperatures, might be useful to reveal the pres- 
ence of a quantum pump effect. 
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APPENDIX 



Eg. 1581 1 to work as a pump is: 



Ring with enclosed magnetic flux: Analytical 
expressions 



The stationary scattering matrix So for the ring with 
branches of length L\ and L 2 , and with enclosed magnetic 
flux $ coupled to two leads via wave splitters S L and S R 
(see, FigtTl reads: 



S = 



So,LL So,LR 
So.RL So,RR 



S ,ll = Sh + St 2 W 2 + St 3 W 3 e L2 (^^ +k ) , 
So.rl = S&W 4 + S? 3 Wie iLl (^^ +k ) , 
So.lr = St 2 Wi + S^W{e L2 (^^ +k ) , 

So,rr = Sg + S? 2 Wi + S&W£e 



nuiA - a 2 StSB , 



&rl &rr 



a M = S&SW* + Sf 3 (sWi -ZTT^ 
a, LR = S&6Wi + Sf 3 (SWI - ^-^) eM 2 ?^ 



di 



Here wz,, ~ v/Li, j = 1,2. The vector-columns 5W 7 and 



5W 



W 1 = -i%M'- 1 %^. 



Note that all the quantities depend on energy E = h 2 k 2 / (2m) 
through the phase factors e L i k , j = 1,2 and on time £ 
through the scattering matrices of wave splitters S^(t) = 
S Wa) (e (t)), (3 = L,R [see, Eqs.gSI]- 



Here L = Li + L2; $0 = ft/e is the single-electron magnetic 
flux quantum. The vector-columns W and W 1 axe defined in 
the following way: 



W = M _1 F, 



W = M'~ X Y\ 



Y 



i 



M = 



M' 



\ 



-1 


V S&c[Li,k] 

( -1 


S£ 3 c[L 2 ,k] 








S32 


-c[Li, - 



oii 

oii 
'-'32 



-c[L 2 ,- 



( S& \ 



Y' 



q R 






/ 



S&c[L2,k] 
Skc[L 2 ,k] 

-1 

k] 

SSc[Li,fc] 
SSc[Li,fc] 

-1 

■k] 



° \ 

-c[L 2 , —k] 



qR 
■^22 





-c[Li, - 

s 22 



) 



Here we introduced the functions c[x, y] = e v L +H / . The 
anti-symmetric matrix A characterizing the ability of the 
ring with oscillating coupling parameters e L (t) and £ R (t) [see, 
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